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We review some previous studies concerning the intra-bilayer Josephson plasmons and present 
new ellipsometric data of the c-axis infrared response of almost optimally doped Bi2Sr2CaCu2 08. 
The c-ajcis conductivity of this compound exhibits the same kind of anomalies as that of underdoped 
YBa2Cu307_5. We analyze these anomalies in detail and show that they can be explained within a 
model involving the intra-bilayer Josephson effect and variations of the electric field inside the unit 
cell. The Josephson coupling energies of different bilayer compounds obtained from the optical data 
are compared with the condensation energies and it is shown that there is a reasonable agreement 
between the values of the two quantities. We argue that the Josephson coupling energy, as determined 
by the frequency of the intra-bilayer Josephson plasmon, represents a reasonable estimate of the 
change of the effective c-axis kinetic energy upon entering the superconducting state. It is further 
explained that this is not the case for the estimate based on the use of the simplest "tight-binding" 
sum rule. We discuss possible interpretations of the remarkable agreement between the Josephson 
coupling energies and the condensation energies. The most plausible interpretation is that the 
interlayer tunneling of the Cooper pairs provides the dominant contribution to the condensation 
energy of the bilayer compounds; in other words that the condensation energy of these compounds 
can be accounted for by the interlayer tunneling theory. We suggest an extension of this theory, which 
may also explain the high values of Tc in the single layer compounds Tl2Ba2Cu06 and HgBa2Cu04, 
and we make several experimentally verifiable predictions. 



I. INTRODUCTION 



The interlayer tunneling (ILT) theory provides a simple explanation of the surprisingly high values of Tc in 

the cuprate superconductors. It is based on the idea ||l|,||j3[^ that the pairing mechanism is substantially amplified by 
a decrease of the effective c-axis kinetic energy of the electrons upon entering the superconducting state. A prerequisite 
for this decrease is the absence (or at least a strong suppression) of the coherent single particle tunnehng between 
the copper-oxygen planes. The onset of the Cooper pair tunneling by the Josephson mechanism at Tc then leads to 
the decrease of the kinetic energy. According to the ILT theory, this gain of energy, which can be expressed as 
the negatively taken coupling energy of the internal Josephson junctions (Ej), represents the dominant part of the 
condensation energy of the superconductor {Uq), 

Ej^Uo. (1) 

This relation, which should be exact in the limit of negligible in-plane contribution to the condensation energy and 
negligible coherent single particle tunneling, has been shown to be only moderately violated for La2_j;Srj;Cu04 [y| but 
strongly violated for TlaBaaCuOg (Tl-2201) [|[0|,§ and for HgBa2Cu04 (Hg-1201) Q. Here we show that Eq. (1) 
is fulfilled for two compounds that have two copper-oxygen planes per unit cell (bilayer compounds): YBa2Cu307_5 
(Y123) and BiaSrzCaCuzOg (Bi-2212). 

So far the relation (1) has been tested only for single layer compounds, since it was not clear what kind of coupling 
takes place for the closely spaced copper-oxygen planes of the bilayer compounds. Anderson assumed that even these 
planes are only weakly coupled and argued that the bilayer compounds consist of two kinds of Josephson junctions: 
inter-bilayer and intra-bilayer | pl| . The Josephson plasma frequency of the inter-bilayer junction {uint) is lower 
than that of the intra-bilayer junction (uJbi). The former determines the c-axis penetration depth while the latter 
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determines Uq. Van der Marel and Tsvetkov proposed a phenomenological model of the dielectric response of 
such a superlattice of inter-bilayer and intra-bilayer Josephson junctions. They showed that it exhibits a transverse 
resonance between the two zero crossings corresponding to the two plasmons. It has been suggested jl^ that this new 
excitation ("transverse plasma excitation"), which can be visualized as a resonant oscillation of the condensate density 
between the two closely spaced copper-oxygen planes, has indeed been observed as an additional absorption peak which 
appears at low temperature in the spectra of the infrared c-axis conductivity of underdoped Y123 []l^ , ^p^|Jl7|] . Very 
recently, this interpretation has been put on a firm basis by a detailed analysis of the c-axis conductivity data for 
Y123 with different oxygen concentrations Note that the observation of the transverse plasma excitation also 

implies the existence of the second (intra-bilayer) Josephson plasmon which is a vital ingredient of the ILT theory. 
This opens a possibility to test also for bilayer compounds the relationship between the Josephson coupling energy 
and the condensation energy predicted by the ILT theory (Eq. (1)). 

In this paper we review some of the previous experimental observations of the transverse plasma excitation (TPE). 
In addition, new data for almost optimally doped Bi-2212 with Tc ~ 91K are reported and analyzed (section II). In 
section III, the Josephson coupling energies of different bilayer compounds obtained from the c-axis conductivity data 
of Refs. |l^jl^,^ , ^ and of the present work are compared with the condensation energies obtained from the specific 
heat data ]22| , |23| . Section IV contains a discussion of the frequently used sum-rules-based estimates of the changes of 
the c-axis kinetic energy |2j,^.^ . Finally, an extension of the interlayer tunneling theory is proposed in section V, 
which allows one to explain the high values of Tc in the single layer compounds Tl-2201 and Hg-1201. 



II. RESULTS AND DISCUSSION 



Since the early days of the high-temperature superconductivity, it was known that some of the infrared active 
c-axis phonon modes of the cuprate superconductors exhibit changes (so-called "phonon anomalies") in the vicinity 
of the superconducting transition temperature (see, e.g. Refs. p7| , p8| , p9| | ) . The most pronounced anomalies have been 
observed for underdoped Y123 0Jll,0,|lll. 



A. YBazCuaOr-a 



The anomalies are illustrated in Fig. 1, which shows our experimental spectra of the c-axis conductivity of 
YBa2Cu3 6.55 with Tc = 53 K jlj]. Note the following anomalies. (A) As the temperature decreases, the oxy- 
gen bond-bending mode at 320 cm~^ which involves the in-phase vibration of the planar oxygens against the Y-ion 
and the chain ions |3l|j3^ , |3^ loses most of its spectral weight and softens by almost 20 cm^^; (B) at the same time 
a new broad absorption peak appears in the spectra around 410 cm~^ and (C) the spectral weight of the peak corre- 
sponding to the apical oxygen mode at 560 cm^^ decreases. The additional absorption peak (feature (B)) has been 
sometimes considered to be a new phonon On the other hand, van der Marel and coworkers suggested that it 
corresponds to the TPE of their model [Q. We believe that this interpretation is correct, but one has to keep in mind 
that some part of the spectral weight of the peak 50%) indeed comes from the phonons (see Ref. for a detailed 
discussion) . With increasing doping the additional peak shifts towards higher frequencies and it becomes broader and 
less pronounced [ p^JT?) . Although the anomalies start to develop above Tc in strongly underdoped compounds, the 
pronounced and steep changes always occur right at Tc Similar though less spectacular effects have also been 

observed for other values of S 16 17, p^ , [l9| , ^ , for several other bilayer compounds | |35| , |36| , ^ , ^ 39|, and for the 

trilayer compound Tl2Ba2Ca2Cu30io [40|. Also shown in Fig. 1 are the fits obtained using the model of Ref. [ p^ . 
Since we are going to refer to this model rather frequently, we summarize its main ideas in the appendix A. It can be 
seen that the model is capable of providing a good fit of both the normal and the superconducting state data. The 
same values of the oscillator strengths of the phonons have been used for the normal and for the superconducting 
states. This means that the changes of the spectral weight are accounted for by the model rather than simply fitted. 



B. Bi2Sr2CaCu208 



Very recently, Zelezny et al. pO|pl[ | have reported similar anomalies in the infrared c-axis conductivity of underdoped 
Bi-2212. The spectra exhibit an increase of the electronic background around 450 cm~^ below Tc accompanied by 
characteristic phonon anomalies. The increase of the background corresponds to the additional absorption peak such 
as observed in the spectra of underdoped Y123. The most pronounced phonon anomaly again consists in a sizeable 
decrease of the spectral weight of the 355 cm~^ phonon mode, which has — according to the shell model calculations 
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PI — a similar eigenvector as the 320 cni~^ mode in Y123. The doping and temperature dependence of the anomahes 
is close to that observed in Y123. These findings of Zelezny et al. confirmed that the anomalies are a common property 
of the bilayer cuprate compounds. 

Fig. 2 shows our experimental spectra of the c-axis optical conductivity of an almost optimally doped Bi-2212 single 
crystal with Tc = 91 K which have been obtained by ellipsometric measurements (see Rcf. [ pT| and references therein 
for a description of the technique). The measurements have been performed on a large ac- face of size 4 x 0.5 mm^. 
As shown in Fig. 2 (b) the spectra exhibit the same kind of anomalies as underdoped Y123 Jl^ , p^ , p^ 17 1^ and 



underdoped Bi-2212 |20| , |2l|| . Below Tc the electronic background increases in the frequency region around 550 cni^^ 
(feature (B)) and simultaneously some of the infrared-active phonon modes are renormalized; in particular, the phonon 
mode at 355 cm~^ loses a large part of its spectral weight (feature (A)). As compared to the data for underdoped 
Bi-2212, the region of the spectral weight increase is shifted towards slightly higher frequencies. This can be expected 
because the frequency of the TPE is known to increase with increasing hole doping (see Ref . ||l^ for a discussion) . The 
onset temperature of the anomalies coincides with Tc (cf. the insets of Fig. 2 (a) and Fig. 2 (b)) which confirms that 
they are related to superconductivity. Note that some signatures of the anomalies can also be found in the reflectivity 
data of Ref. ||. 

Figure 3 (a) displays the region near the anomalies on an enlarged scale and Fig. 3 (b) shows the fits obtained by 
using a slightly modified version j2^ of the model of Ref. ||l^ . The modification concerns only the description of the 
featureless electronic background. First, the inter-bilayer susceptibility is set equal to zero {xint = in notation of 
Ref. ) . This is a reasonable approximation considering the very low value of the unscreeneed inter-bilayer plasma 
frequency Uint {^int < 20cm~^, see Ref. [^). Second, the regular part of the intra-bilayer susceptibility [xu) bas 
been fitted by a combination of a broad Drude term and a broad mid-infrared Lorentz oscillator instead of the single 
very broad Lorentzian centered at high frequencies which was used previously. Furthermore, we have adopted the 
simple picture of the phonon eigenvectors used in the fits of the spectra of underdoped Y123 and underdoped Bi-2212, 
i.e., we assume that the 355 cm~^ mode corresponds to the vibrations of the planar oxygens, whereas the modes at 
300 cm~^ and 580 cm~^ correspond to vibrations of the inter-bilayer oxygens, i.e., the apical oxygens and the oxygens 
of the BiO layers. The values of the fitting parameters and other numerical factors used are summarized in Table 
I. Those corresponding to the room-temperature spectrum have been obtained by fitting the measured frequency 
dependence of the complex dielectric function in the frequency region above 300 cm^^ (with lou = and Su =0, Su 
being the oscillator strength of the mid-infrared Lorentz oscillator). Those used in calculating the 100 K spectrum 
have also been obtained by fitting the data (with lou = and Su = 0), except for the values of Eoq, 5'p, 5*1, ^2 
(oscillator strengths of the phonons), ojp (frequency of the oxygen-bond bending mode), lo\ and 71 (parameters of 
the 300 cm^"'^ mode), which have been fixed at the room-temperature ones. Finally, the values of £00, S'p, 5*1, 5*2, 
^Pi 7Pi ^1i71i E^nd (broadening parameter of the broad Drude contribution), have been fixed at the lOOK ones, 
when fitting the low-temperature spectrum. It can be seen that the most pronounced features of the experimental 
spectra, in particular the increase of the electronic background around 550 cm^^ and the spectral- weight anomaly of 
the oxygen bond-bending mode, are well reproduced. 

Figure 3 (c) shows the frequency dependence of the difference CTi(r = lOK << Tc) — a\{T = 100 K w Tc). The 
value of the difference is positive in the frequency region between 420 cm~^ and 580 cm~^, whereas it is negative both 
for lower frequencies (in the region between 340 cm~^ and 420 cm^^) and for higher frequencies (in the region between 
580 cm~^ and 650 cm~^). The positive values in the first region are caused mainly by the increase of the electronic 
background below Tc (i.e., the additional peak due to the TPE). The negative values in the second region correspond 
to the spectral weight anomaly of the oxygen bond-bending mode. The negative values in the third region are mainly 
due to a slight shift of the apical-oxygen mode at 580 cm~^ towards lower frequencies. The agreement between the 
difference of the measured data (full line) and the difference of the fitted spectra (dashed line) is excellent. Note, that 
this agreement has not been achieved by changing the values of the parameters of the phonon modes. The slight shift 
of the 580 cm~^ mode accounts only for a part of the spectral weight increase around 550 cm~^. 

Finally, the frequency dependence of the quantity Nn — Ng, where Nn{Lj) — N{T = 100 K fa Tc,uj), 
Ns{uj) = N{T = lOK « Tc,w) and 



iV(T,w)= / duj'ai{T,uj'), (2) 
Jo+ 

is displayed in Fig. 3 (d). We take our low-frequency cutoff of 70 cm"'^ as the lower limit of the integration in Eq. (2), 
but it is rather unlikely that below this frequency there is any considerable difference between the normal state and 
the superconducting state data. For a conventional superconductor, the value of A^^^ — Ng increases with increasing 
frequency and it approaches ps, the spectral weight of the 5- peak at w = 0, within a range of ~ 6 A, A being the 



superconducting gap |45|. It has been shown by Basov et al. [ |25[| that for several high-temperature superconductors 



the value of Nn — Ng also increases with increasing frequency but it saturates at a value of only about one half of 
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Ps for Lo ^ lOA. Our result for Bi-2212 exhibits an even more suprising tendency: Nn — Ng stays approximately 
constant up to ~ 350 cm~^, it changes considerably in the frequency region between 350 cm~^ and 620 cm~^ and it 
seems to saturate above ~ 650 cm~^ at a negative value of approximately — 300 f2^^ cm^^. The corresponding value 
of the ratio |A'^,i — Ns\/ps is larger than 25. The fact that the value of Nn — Ng above ~ 550 cm^^ is negative signals 
an increase of the low frequency spectral weight below Tc which is caused by the formation of the peak due to the 
TPE. Note that such an increase cannot be easily explained by any conventional theory. We shall come back to this 
interesting issue in section IV B. 



III. JOSEPHSON COUPLING ENERGIES AND CONDENSATION ENERGIES 



In Table II we summarize the values of the intra-bilayer plasma frequency (ujbi), the corresponding Josephson 
coupling energy per unit cell (Ej), and the condensation energy per unit cell {Uq) for Y123 and Bi-2212. The values 
of Ej have been calculated using the formula p6| 

Ej = \ or equivalently £;j[meV] = — ^j- (wfc;[cm"^])^ . (3) 

Here a is the in-plane lattice constant, du is the distance between the closely spaced copper-oxygen planes, and 
C = 3.1 X 10^^. We neglect the contribution of the inter-bilayer plasmon to Ej which is by at least one order of 
magnitude lower than that of the intra-bilayer one. Note that Eq. (3) does not contain any adjustable parameters. 
The values of the condensation energies have been obtained from the specific heat data in Ref. (Y123) and in 
Ref. (Bi-2212). We estimate the error bars of the values of uju to be about 20 — 30%. These error bars arise 



mainly from the uncertainty in the description of the electronic background |47|. The corresponding uncertainties of 
the values of Ej are about 50%. Taking this into account, the agreement between the Josephson coupling energies 
obtained from the optical data and the condensation energies obtained from the specific heat data is reasonably good: 

(i) The values of the two quantities are of the same order of magnitude. 

(ii) Both quantities exhibit a rather similar dependence on doping. It is even possible to understand, why the values 
of the ratio Ej /Uq are higher for the strongly underdoped samples than for the less underdoped ones or the optimally 
doped ones. It is namely likely that for strongly underdoped samples some fluctuation effects set in well above the 
macroscopic transition temperature T^. This suggestion is consistent with the finding that the phonon anomalies 
start to occur well above in strongly underdoped samples ||l^,|l^,Q . The contribution of the fluctuation effects to 
the condensation energy is not likely to be contained in the values of Uq presented in Table II, since they have been 
obtained by an analysis of the specific heat data where only the changes occuring below are properly taken into 
account. On the other hand, the values of Ej do contain the contribution of these effects, since they are determined 
by the low-temperature values of uju (i.e., in a way which does not require any assumptions concerning the onset of 
superconductivity). As a result the values of C/q for strongly underdoped samples can be expected to be smaller than 
the values of Ej. 

(iii) The values of Ej for Y123 are systematically higher than those for Bi-2212, in agreement with the trend in the 
condensation energies. It may be argued that the difference in the condensation energies is related to the presence of 
metallic chains in Y123. However, the chain condensation cannot fully account for the difference, since Ca substituted 
samples with broken chains also have a significantly higher condensation energy than Bi-2212 p^ . 

This agreement means that the values of the condensation energies of the bilayer compounds can be explained 
using the ILT theory. Further implications will be discussed in Section V. It certainly remains an open question why 
different compounds with almost identical values of Tc have rather different values of the condensation energy. 



IV. SUM RULES 



There is another approach to estimate the changes of the c-axis kinetic energy, pioneered by Chakravarty, Kee 
and Abrahams which is based on the use of the optical sum rule. For the bilayer compounds, this approach 

yields an opposite result, namely that the value of the condensation energy cannot be accounted for by the change 
of the c-axis kinetic energy. In the first part of this section we clarify several conceptual points which have made the 
discussions of the kinetic energy changes somewhat confusing (subsection A) . Next we discuss some specific properties 
of the bilayer compounds (subsection B). In particular, we show that for these compounds the change of the c-axis 
kinetic energy cannot be estimated using the conventional "tight-binding" sum rule of Refs. 



24|. 
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A. Conceptual issues 



For simplicity, we focus in our considerations first on the single layer compounds. In the present context, there 
are at least three different quantities denoted as c-axis kinetic energy. First, the "true" c-axis kinetic energy, Kc — 
^iPi.cf^i^- Second, the "tight-binding" c-axis kinetic energy He, 

Hc^~J2 t^il^hj)4j+,,,c,,Ls + H.c. , (4) 

l,i,j,s 

where 1 is the layer index, i and j refer to the sites of the two-dimensional layers, s refers to spin, t± is the interlayer 
hopping matrix element, and H.c. means the Hermitian conjugate operator. Third, the effective low-energy c-axis 
kinetic energy of the ILT theory Hj, 

Hj ^ - Yl rj(k)Ck,i+i,tcik,i+i,i'=-k:i,iCk,i,T + H.c. , (5) 

i.e., the Josephson coupling energy. The experimental data are discussed in terms of at least two different sum rules. 
First, the "general" sum rule (for derivation see, e.g, Ref. ||49|]), 

/■°° / ^ ^ne'^ 

which is fulfilled for any real system and second, the "tight-binding" sum rule valid only for model Hamiltonians 
whose single-particle part is of the tight-binding form (for derivation see, e.g., Refs. ||5^,^ and references therein), 

r a,{T,u;)<h;^^^{He)iT). (7) 

In Eq. (6), n is the total electron density, in Eq. (7) a is the in plane lattice constant, d is the lattice constant along the 
c-axis, and He is the tight-binding kinetic energy (Eq. (4)) per unit cell {{He) represents an average of this quantity, 
T is the temperature). 

Let us discuss possible changes of the three kinetic energies upon entering the superconducting state and the 
relations between them and the sum rules. 

(i) There is no obvious relation between {Kc) and the sum rules. Upon entering the superconducting state, the total 
energy of the superconductor decreases. This implies — by virtue of the virial theorem {{K) = —{V)/2, where K and 
V are the kinetic and the potential energy, respectively; see, e.g., Ref. |^l[]) — that its total kinetic energy increases. 
It is possible to assume that even {Kc) increases at the superconducting transition. 

(ii) It is very likely that the degrees of freedom essential for superconductivity are contained in some model single 
band Hamiltonian Hu, whose (tight-binding like) single particle part is derived by a downfolding process in which 
all of the higher energy bands are integrated out ||5^ and whose interaction part describes the most pronounced 
correlation effects. In our opinion, there are no general reasons why the average of its c-axis kinetic energy term He 
(Eq. (4)) should decrease or increase below Tc. The arguments invoking a frustration of the c-axis kinetic energy 
in the normal state and a "deconfinement" of bosonic pairs below Tc (See, e.g., Ref. Chapter 2, Dogma VI and 
Ref. IH) are usually based on asymptotic properties of the quasiparticle propagators, i.e., they concern rather the 
asymptotic behaviour for lu —^ than the actual changes of {He). Possibly these arguments could be formulated more 
precisely in terms of yet another effective kinetic energy H'c characterising the c-axis electrodynamics in a frequency 
scale intermediate between that of He (~ 10000 cm~^) and the frequency scale of Hj of the order of 100 cm~^ |5^] 
(see Fig. 4). 

It follows from Eq. (7) that the possible change of {He) at the superconducting transition is related to the change 
of the spectral weight as follows: 

{{He)s{T «Tc)- {He)n{T ~ T,)) [ mcV] = 



4C 120 / 



d[k] TT 



ai{T «Tc,uj)duj- (Ji{T KTc,uj)duj\ [n-^cxar'^], (8) 



where VLc is a cutoff frequency required to exhaust the sum rule (7). Note that in a model involving strong correlations, 
the value of fig may be much higher than the bandwidth of the corresponding noninteracting model. In addition, this 
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value (and also the changes of the kinetic energy themselves) may depend on the kind of model which is used, i.e., on 
the interaction part of the model Hamiltonian Hu,. 

(iii) We emphasize that it is the effective c-axis kinetic energy Hj of the ILT theory (Eq. (5)) whose changes at the 
superconducting transition were predicted to be responsible for the high values of Tc in Ref. . The only "sum rule" , 
which can be associated with Hj, is Eq. (3) with du substituted by d (note that Ej = —{Hj)). 

What is the experimental status? The analysis of Basov and coworkers |^ demonstrates that at least for the 
underdoped high- Tc superconductors, the low- frequency spectral weight 

a{T,cu)= / ai{T,u;')du;' ^p,{T)+N{T,Lj) (9) 

increases considerably upon entering the superconducting state for cutoff frequencies oj ranging up to at least 
1200 cm~^. More precisely, as one increases the value of uj, the value of the difference a{T << Tc,uj) — a{T « Tc, 
saturates at « 500 cm~^, reaching typically 50% of ps- The only way to reconcile this interesting finding with the 
general sum rule (6) consists in assuming [g5[ that a substantial part of a{T « Tc,1200cm~^) is collected from 
frequencies exceeding 1200 cm~^. In other words, since the total spectral weight has to be conserved, the increase 
of a(r, 1200 cm~^) below Tc is compensated by a decrease of spectral weight around some frequency lu* higher than 
1200 cm^^. A sketch of the situation is shown in Fig. 4. 

What do we learn from these results about the changes of the tight-binding kinetic energy {Hc)l Let us denote by 
LOc the interband cutoff, i.e., the upper limit of the frequency interval for which the response of the superconductor 
can be properly described using a certain effective Hamiltonian Hu, (e.g., the Hubbard Hamiltonian). Let us further 
assume, for the sake of simplicity, that lOc > of Eq. (8) There are two different possibilities, which are both 
compatible with the infrared data: (a) lo* > ujc and (b) uj* < ujc (see Fig. 4). In case (a), it follows from Eqs. (6) 
and (8) that (He) indeed decreases below Tc- In order to establish this decrease experimentally, however, one would 
have to estimate the total low-frequency spectral weight above and below Tc, i.e., to integrate ai{uj) up to ujc- We 
believe that the case (b) is more likely to happen. The increase of a{T, 1200 cm~^) below Tc and the corresponding 
decrease of the spectral weight around uj* then both can be described by using Hu,. It follows from Eqs. (6) and 
(8) that (He) does not change at Tc- To conclude, it is only the tight- binding c-axis kinetic energy He (Eq. (4)) of 
an effective single band Hamiltonian, whose change at the superconducting transition can be obtained using Eq. (8). 
The results of Basov et al. nicely reveal the unconventional properties of the high-temperature superconductors, 
in particular the extremely large frequency scale involved, but they cannot be used to yield reliable estimates of this 
change because of the uncertainties concerning the cutoff frequencies flc and cUc- In addition, there are no apriori 
reasons, why (He) should change upon entering the superconducting state. 

Another example of a model c-axis kinetic energy is represented by the effective low-energy c-axis kinetic energy 
Hj (Eq. (5)) of Ref. 0. It acquires a nonzero (negative) value, —Ej (Eq. (3)), only in the superconducting state and 
it was predicted to be responsible for the high values of Tc- It is certainly a crude approximation to identify Ej with 
the c-axis contribution to the condensation energy, as we did it in the previous section, since possible "countereffects" 
p5| are not included. They are to some extent included when estimating the changes of the c-axis kinetic energy 
using Eq. (8) with the upper limit of the integrals well below lu* (instead of Wc), i.e., 

AC 120 

AEk^nA^eV]^^—J- {a{T «Tc,uj<u;*)-a{T^Tc,iu<u;*)) [fl^^ cir-^] , (10) 

d[A\ TT 

This is an appealing possibility, but the physical meaning of the result is not completely clear (in contrast to the 
change of the well defined quantity (He))- It may be perhaps viewed as an estimate of the change of some kinetic 
energy H'c characterising the c-axis electrodynamics in an intermediate frequency scale (see our discussion above 
and Fig. 4). Note finally that the fact that AEkin,c includes the countereffects does not necessarily imply that 
\^Ekin c\ < Ej because of the additional factor of 4 entering the sum-rules-based formulas [Q. On the contrary, for 
moderate countereffects, \AEkin,c\ > Ej. For underdoped La2-xSr2,Cu04, e.g., \AEkin,c\ is by a factor of 2 larger 
than Ej E^. 



B. Some specific properties of the bilayer compounds 

In Fig. 3 (d), we have encountered a rather unusual situation. First, the value of the quantity A'^(r, 690 cm^^) 
defined by Eq. (2) increases with decreasing temperature. A phenomenon not observed for any of the single-layer 
compounds (see, e.g., Ref. [|5|). Second, the increase of the integrated spectral weight a(r, 690cm~^) defined by 
Eq. (9) upon entering the superconducting state is lower by a factor of ~ 30 than the value required to yield a kinetic 



6 



energy change comparable to Uq (when inserted into Eq. (10)). On the other hand, using simply the formula (3) for 
the Josephson coupling energy, we have obtained a value fairly close to Uq. Here we propose a qualitative explanation 
of these paradoxes. 

Let us consider the superlattice of intra-bilayer and intcr-bilayer Josephson junctions. The low frequency spectral 
weight a{T << Tc,uj > ujp), where ujp is the frequency of the TPE (Eq. (A. 2)), contains both the contribution of the 
i5-peak at a; = 0, Ss (Eq. (A. 3)), and the contribution of the TPE at a; = ujp, Spi (Eq. (A. 4)). The appearance of the 
second one below Tc explains the first paradox. It is easy to show that 

a{l « lc,uj > Up) - -Sq , , , , — , (llj 

where dint is the distance between the bilayers. The Josephson coupling energy is approximately given by Eq. (3) 
(neglecting the contribution of the inter-bilayer Josephson junction), whereas the change of the c-axis kinetic energy 
estimated using Eqs. (10) and (11) is given as follows: 

^J. r 4C du (a;M[cm~i])^ 

A£;fc„,c[meV] = -— — — — ' , , — ■ (12) 

[dbl + dint ) [AJ dbl + dint 

It can be seen that 



EJ»\AEt^n,c\. (13) 

This is the explanation of the second paradox. The effects of the electronic background make the discrepancy even 
more pronounced. At the same time, it is obvious that within the simple model of the superlattice of intra-bilayer 
and inter-bilayer Josephson junctions, it is the Josephson coupling energy Ej, rather than AEkin,c, which represents 
the change of the c-axis kinetic energy upon entering the superconducting state. We conclude that for the bilayer 
compounds the change of the c-axis kinetic energy can be estimated using Eq. (3) but cannot be estimated using 
Eq. (10). We refer the reader to appendix B for a more rigorous discussion. 



V. POSSIBLE EXTENSION OF THE ILT SCENARIO FOR SINGLE LAYER COMPOUNDS 



In section III we have shown that for the bilayer cuprate compounds the Josephson-coupling energy such as estimated 
from the optical data can account for the condensation energy. In section IV B we have shown that this result is not 
invalidated by the fact that Eq. (10) yields a much lower estimate of the change of the c-axis kinetic energy upon 
entering the superconducting state than Eq. (3). The main reason of the discrepancy is that the tight-binding sum 
rule of Eq. (7) used to derive Eqs. (8) and (10) is only valid for such single layer compounds, for which the distribution 
of the total electric field is homogeneous. The results presented in Section III thus suggest that the high-tcmpcrature 
superconductivity in the bilayer compounds can be accounted for by the ILT theory. We are left with three possible 
explanations. 

(a) The agreement between the Josephson-coupling energies and the condensation energies reported in section III 
represents a mere coincidence. 

(b) The interlayer tunnelling indeed provides the dominant contribution to the condensation energy of the bilayer 
compounds. Another mechanism is responsible for the high values of Tc in the single layer compounds Tl-2201 and 
Hg-1201. 

(c) A modified ILT theory may explain the high-temperature superconductivity both in the bilayer compounds and 
in the single layer compounds. 

In our opinion, it is unlikely that case (a) is realized. It would mean either that the assignment of the additional 
absorption peak to the TPE |l^,|l^,|l^ is wrong or that the Josephson-coupling energy is completely (or almost 
completely) compensated by the countereffects There are several important arguments supporting the present 
interpretation of the data. Let us mention two of them. First, the doping dependence. The frequency of the 
maximum of the additional absorption peak increases considerably with increasing doping. This is easy to explain 
within the present scenario (the squared frequency of the maximum should be proportional to the condensate density) 
but difficult to explain within theories where the additional peak is assigned to a pair breaking excitation p7| , ^ . 
According to these theories the frequency of the maximum should be close to the pair breaking frequency observed 
in other experiments, which seems not to be the case for strongly underdoped samples. One could argue |59|] that the 
final state interactions may shift the peak towards lower frequencies. Second, both the frequency and the spectral 
weight of the additional peak are determined by a single parameter, oju ||60| . In other theories, at least two parameters 
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are required to fit the data, e.g., A and t±. Starting from our interpretation of the data, large countereffects would 
represent the only possibility to rule out the ILT mechanism from the role of the mechanism providing the dominant 
contribution to the condensation energy Uq. Note that the values of the Josephson-coupling energies presented in 
Table II are so high that even sizeable countereffects of the order of 80 % would not invalidate our interpretation (see 
Ref. [^). A possible analysis of the countereffects along the lines of appendix B is complicated by the fact that the 
normal-state data for very low temperatures, which should be compared to the superconducting-state data, rather 
than the data obtained for temperatures above Tc, are not available. If case (b) were realized, there would be two 
different mechanisms causing the high values of Tc- This is certainly not impossible. In the remaining part of this 
section we argue that even the possibility (c) should not be excluded. We propose an extension of the ILT theory, 
which may explain the high- Tc values in the single layer compounds Tl-2201 and Hg-1201. For concretness we focus 
on the Tl-2201 compound, for which some optical data are already available P,P,p|, |6l| , p5p^ ] . 

According to the ILT theory as formulated, e.g., in Ref. the c-axis kinetic energy is frustrated in the normal state 
because of the spin-charge separation mechanism. The electrons or holes are composite objects that cannot escape from 
the copper-oxygen planes (so called "confinement" ) . For Cooper pairs this confinement is relaxed and consequently 
the c-axis kinetic energy decreases at the superconducting transition. To our best knowledge, it has been previously 
always assumed that this decrease is related to the onset of Josephson tunneling between the superconducting copper- 
oxygen planes. This is, however, only one possible mechanism for the decrease. We suggest that in Tl-2201 the c-axis 
kinetic energy decreases via some delocalization of the Cooper pairs which involves the apical-oxygen orbitals. By the 
delocalization we mean the fact that the superconducting-state wave function acquires a larger contribution of the 
apical-oxygen orbitals than the normal-state wave function. Quite generally, the c-axis kinetic energy decreases at the 
superconducting transition, whenever the Cooper pairs are more delocalized along the c-axis than the single electrons 
or holes in the normal state. Note that our extension of the ILT theory is based on a radical form of the "confinement" 
hypothesis. Not only the hopping between the copper-oxygen planes is assumed to be blocked in the normal state but 
also (at least to some extent) the hopping between the orbitals of the copper-oxygen planes and the apical-oxygen 
orbitals. The conjecture can be experimentally tested in several ways. First we discuss the c-axis infrared spectra, 
which may already provide some evidence in favour of the proposed scenario. Next we briefly mention some other 
experimentally verifiable consequences. 

Recently, some infrared data for Tl-2201 have been reported [^,^. The spectra of the c-axis conductivity exhibit 
four peaks corresponding to the c-axis infrared phonons — instead of five as dictated by symmetry — and two other 
distinct features, which have not been discussed previously: 

(i) A pronounced step around 210 cm^^. For convenience, the relevant part of the conductivity spectra is shown in 
the inset of Fig. 5 (b). Note that a weak feature at ~ 230 cm"-'^ in the grazing incidence reflectivity data of Ref. |6|] 
may be related just to this step. We suggest that the step is due to the missing fifth c-axis infrared phonon. 

(ii) An increase of the spectral weight of the phonon peak at 390 cm~^ upon entering the superconducting state. It 
can be concluded from Fig. 3 of Ref. that the spectral weight in the frequency region between 280 cm"-'^ and 
440 cm~^ slightly (but not negligibly) increases below Tc- For convenience, the important part of the data is shown 
in the inset of Fig. 5 (c). We estimate that the spectral weight increase represents 3-5% of the spectral weight of the 
phonon. 

Below we show that both features (i) and (ii) are compatible with the suggested extension of the ILT theory. 

Figure 6 shows the crystal structure of Tl-2201 such as reported in Ref. [Q. The simplest way to explore the 
possible consequences of our hypothesis consists in dividing the structure into regions A and B (denoted in Fig. 6), 
assuming that they exhibit different properties (region B being rather insulating and region A being more metallic), 
and applying the model of Ref. |l^ ]. In order to be able to apply the model, we have to make some assumptions 
concerning the eigenvectors of the phonons. Similarly as in case of Bi-2212 we concentrate on the three phonons 
located at higher frequencies (See Fig. 2 of Ref. [^) presumably corresponding to vibrations of the oxygens: the 
phonon mode at 600 cm~^, the mode at 390 cm~^, and the phonon which shows up as the step at 210 cm~^. According 
to the shell model calculations of Ref. |Q, the highest frequency mode involves vibrations of the TlO-layer oxygens, 
the second one mainly vibrations of the apical oxygens, and the third one vibrations of the planar oxygens. Besides 
this assignment (assignment Al), we have also considered the possibility that the highest frequency mode corresponds 
to vibrations of the apical oxygens (as in many other cuprates) and the second one to vibrations of the TlO-layer 
oxygens (assignment A2). In reality the vibrations of the apicals and the TlO-layer oxygens are probably strongly 
coupled. Figure 5 shows the results of our simulations. Parts (a) and (b) have been obtained using the assignment 
Al, parts (c) and (d) using the assignment A2. The model in the form presented in Ref. allows us to treat only 
two phonon modes at the same time (the "interface" one, i.e., the apical-oxygen mode, and only one of the other 
two modes). This is the reason why there are two figures for a given assignment, each showing only two structures 
related to the phonons. The electronic conductivities of the two regions A and B [ga and as) have been modelled, 
for the sake of simplicity, by broad Drude terms, the plasma frequency of the region A, Q.a, being much higher than 
the plasma frequency of the region B. The results are shown for two different values of VIa'- ^ai (dashed line) and 
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^A2 (solid line), Q,ai < ^A2- The increase of Qa may simulate the changes brought about by the superconducting 
transition: the more delocalized the ground state wavefunction, the more metallic the region A and the higher the 
value of il-A- For the sake of simplicity, we do not consider the changes of as caused by the onset of superconductivity. 
The values of the parameters used are given in Table III. 

It can be seen in Figs. 5 (b) and 5 (d) that the lowest frequency mode manifests itself as a step rather than as 
a Lorentzian peak, in agreement with the experimental data. Note that the step is getting less pronounced with 
increasing value of Qa, i.e., with increasing metallicity of the region A. This may correlate with an apparent absence 
of the structure in the data for strongly overdoped Tl-2201 (See Fig. 2 of Ref. fG^])- As seen in Figs. 5 (a) and (c), 
the phonon mode involving vibrations of the TlO-layer oxygens exhibits distinct changes with increasing value of O^, 
in particular its spectral weight increases |]65|| . This may correspond to the observed spectral weight anomaly of the 
390 cm^^ phonon mode (assuming that the assignment A2 is closer to reality than the assignment Al). Finally, the 
increase of JIa results in an increase of the electronic background around 1000 cm"-'^ (concerning the absolute values 
see Ref. [p6|). Such an increase should be observed in future experiments. Note that in the present simulations, 
the increase of ^Ia is intentionally relatively high so that the anomalies are clearly seen. Such an increase would 
correspond to a change of the c-axis kinetic energy much higher (by a factor of ~ 4, a rough estimate obtained by 
using Eq. (B.ll)) than the condensation energy |67| ]. The actual changes of the electronic background to be observed 
experimentally therefore should be considerably smaller (e.g., by a factor of 4). Also the frequency range of the 
increase may be somewhat different. To summarize, our simulations demonstrate that the anomalous features (i) and 
(ii) can be explained within a model involving the variations of the electric field inside the unit cell, whereas they 
cannot be easily explained in a conventional way. The important point is that if the electric field indeed changes 
considerably inside the unit cell, the changes of the c-axis kinetic energy upon entering the superconducting state 
cannot be estimated simply from the value of the plasma frequency of the superconducting condensate or by using 
the modified tight-binding sum rule (10). 

Let us compare the c-axis conductivity of Tl-2201 with that of another single layer compound with considerably 
lower Tc, Bi2Sr2Cu06 (Bi-2201). The spectra for Bi-2201 exhibit five peaks corresponding to the c-axis infrared 
phonons (See Fig. 2 of Ref. |6|], similar results have been obtained by our group ^^), i.e., there is one phonon peak 
more than in the spectra for Tl-2201. This difference finds a natural explanation within our extension of the ILT 
theory. Bi-2201 has a substantially lower condensation energy than Tl-2201 ||6^,|7^. This suggests that in Bi-22G1 
the apical-oxygen orbitals are less accessible for the charge carriers. Consequently, the region A is much less metallic 
than in case of Tl-2201 and the planar-oxygen mode manifests itself rather as a peak than as a step feature. This is 
demonstrated in Fig. 5 (d), where the results for ^Ia = 1500 cm" ^ (dotted line) and Ha — 600 cm" ^ (dashed-dotted 
line) are shown. 

An observation of the infrared anomalies discussed above in more precise future experiments would represent a 
clear but indirect evidence in favour of the suggested extension of the ILT theory. Direct evidence could be provided 
by experiments probing the occupation of individual atomic orbitals, as, e.g., the near x-ray absorption fine structure 
(NEXAFS) measurements (see, e.g., Ref. [0). We predict that such experiments will reveal a sizeable increase of 
the number of holes on the apical sites at Tc accompanied by the corresponding decrease of the number of holes in 
the copper-oxygen planes. Note that this prediction does not concern only Tl-2201 but all the single layer high- Tc 
cuprates. It may, of course, also be fulfilled for some of the bilayer cuprates. The changes of the hole distribution 
at Tc should further result in some structural changes, in particular, in changes of the distance between the apical 
oxygens and the copper-oxygen planes, and in anomalies of some Raman-active phonon modes involving vibrations 
of the apical oxygens. Such anomalies may already have been observed [Q. 

VI. SUMMARY 

The c-axis infrared conductivity of optimally doped Bi-2212 exhibits the same kind of anomalies as that of un- 
derdoped Y123. Below Tc the electronic background increases in the frequency region around 550 cm"^ and at the 
same time the oxygen bond-bending mode at 355 cm~^ loses a large part of its spectral weight. The anomalies can be 
explained within a model involving the intra-bilayer Josephson effect and variations of the electric field inside the unit 
cell. We have compared the Josephson-coupling energies of Y123 and Bi-2212 with different oxygen concentrations 
obtained from the optical data with the condensation energies obtained from the specific heat data and we have found 
that there is a remarkable agreement between the values of the two quantities. The Josephson coupling energy is 
shown to represent a reasonable estimate of the change of the c-axis kinetic energy upon entering the superconducting 
state and it is also shown that the latter quantity cannot be obtained by using the simple "tight-binding" sum rule, 
since this sum rule has been derived assuming homogeneous distribution of the total electric field within the unit cell. 
The most plausible interpretation of the agreement between the Josephson coupling energies and the condensation 
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energies is that the condensation energy of the bilayer compounds can be accounted for by the interlayer tunnehng 
theory. We propose a modification of this theory which may also explain the high values of Tc in Tl-2201 and Hg-1201. 
The main idea is that the c-axis kinetic energy of these compounds decreases at Tc via a delocalization of the Cooper 
pairs onto the apical-oxygen orbitals. We have investigated, using a toy model, the consequences of this hypothesis 
for the c-axis infrared response and we have demonstrated that it offers a simple explanation of two features observed 
in the measured c-axis conductivity of Tl-2201. In addition, we propose that for Tl-2201 an increase of the electronic 
backround around 1000 cm~^ takes place at Tc- We further predict a sizeable increase at Tc of the number of holes 
on the apical sites, related structural changes and related anomalies of some Raman-active phonon modes. 

Note added. There are some similarities between our model |l^ of the charge dynamics and the models of 
Refs. 75 T^JT^. In Rcf. |Q it is argued that the role of the intra-bilayer plasmon could be played by an 

antiphase oscil ation between the copper-oxygen planes and the ordered blocking layer plaques. After submitting the 
paper, O. K. Andersen has drawn our attention to a preprint of Pavarini et al. [7^ , where it is shown that there is a 
correlation between the maximum value of Tc for a given compound, Tcmax, and the value of the parameter r which 
expresses the range of the intra-layer hopping, t' /t, obtained from band structure calculations. The higher the value 
of r, the higher the value of Tcmax- It remains an open question to what extent this correlation is consistent with 
the extension of the ILT theory proposed in Sec. V. There we have suggested that the condensation energy of the 
single layer compounds is largely due to a change of the c-axis kinetic energy at Tc connected with an increase of the 
number of holes on the apical sites. The reason for this change is the confinement of the holes to the copper oxygen 
planes in the normal state and a deconfincment of the Cooper pairs below Tc- In case of the single layer systems 
the condensation energy — and also Tc — thus should be the largest (ignoring the material dependence of the orbital 
energies) for the compounds with the largest value of the hopping matrix element tpa between the relevant orbital of 
a copper-oxygen plane (predominantly dx2-y2) and the 2pz orbital of the neighbouring apical oxygen plane. The 2pz 
orbital of an apical oxygen couples to the copper 4s orbital and tpa is thus determined (a) by the matrix element tsc 
between the 2pz orbital and the Cu 4s orbital and (b) by the ratio of the Cu 4s to Cu dr^2_y2 character, R. The value 
of tsc decreases |78| with increasing distance d between the copper-oxygen plane and the apical-oxygen plane and it 
is thus somewhat smaller for Tl-2201 (d « 2.1 X) and Hg-1201 {d « 2.8l) as compared to LaSrCu04 {d « 2.4A). On 
the other hand, the value of R increases with increasing value of r as fTq], which means that it is much larger for 
Tl-2201 (r « 0.33) and Hg-1201 (r « 0.33) than for LaSrCu04 (r « 0.17) ||7|]. For this reason it is possible that tpa 
increases with increasing r, i.e., that the correlation reported in Ref. [t^ is consistent with what one would expect on 
the basis of the extended ILT theory. Obviously, a careful analysis involving the details of the Cu 4s — apical-oxygen 
2pz bonding is required. 
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Appendix A: The main ideas of the model of Ref. |18| 



Figure. 7 (a) shows a schematic representation of the electronic part of the model of Ref. ||l8|] . The effects related to 
the ionic degrees of freedom will be discussed subsequently. Thick horizontal lines correspond to the two-dimensional 
copper-oxygen planes. The distance between the closely spaced copper-oxygen planes forming a bilayer is denoted by 
du, the distance between the bilayers is denoted by dint- If an electric field E' is aplied, the currents ju (intra-bilayer 
current) and jint (inter-bilayer current) flow between the planes. Since they are not equal, the planes become charged. 
The resulting surface charge density which alternates from one plane to the other is denoted by k. It modifies the 
average electric field in the intra-bilayer region, Eu = E' + (k/eoEoo), whereas it does not change the average electric 
field in the inter-bilayer region, Eint = E'. The currents ju and jint are determined by the fields Eu and Eint, 
respectively, jbi — cruEu and jint — <^intEint, so that a selfconsistent set of equations is obtained. Here (Jbi and 
aint are the intra- and the inter-bilayer conductivities, respectively. The equations can be readily solved yielding the 
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following formula for the macroscopic dielectric function: 



e{u) = {dbl + dint)/ 



+ 



Ebliuj) £int{(^) 



(Al) 



Here £6i(w) = {i/eoUj)abi{Lo) and £mt((jj) = (i/eo'^)o'int(w). For a superlattice of intra-bilayer and inter-bilayer 
Josephson junctions we have eu{(^) — £00 — {'^ul^'^) ^'^'^ £mt(w) = £00 — {'^fnt I ^'^) ■ The dielectric function of 
Eq. (A.l) then exhibits two conventional plasma resonances at w = ^ul \/£ 00 and at cj = <^mt/^£oo- In addition it 
exhibits a pole at 



_ / dbl^lnt + d^nt^ll . 

uj-ujp-i— — — — , {A.2) 

y (dbl + dmt )£oo 

corresponding to a resonant oscillation of the condensate density between the two closely spaced copper-oxygen planes 
(transverse plasma excitation). The spectral weight of the (5-peak at a; = is 



Q {dbl + dint)uJu^int ( a n\ 

dblUJtnt + d^ntUJbl 



and the spectral weight of the resonance at LOp is 



dbldint i^bl 
' {dbl + dint){dblUjfnt + dtnti^b 



Equations (A.l) and (A. 2) have been previously derived in another way by Van der Marel and Tsvetkov (see 
also Ref. |1^). As far as only the electronic degrees of freedom are concerned, the two approaches yield completely 
equivalent results. However, the present approach can be more easily extended to incorporate the phonons. 

Assume for a moment that the displacements of the ions do not modify the electric fields Ebi and Eint ■ The ions 
located in the intra-bilayer and in the inter-bilayer regions experience then the electric fields Ebi and Eint , respectively. 
The ions located in the copper oxygen planes experience a field EiocP which is equal to the average of the two fields, 
EiocP = {Ebi + Eint)/'2. This can be easily explained using Fig. 7 (a). Let us consider, e.g., the ions in the middle 
copper-oxygen plane. They experience the applied field E' and the electric field generated by the charge density —n 
of the upper copper-oxygen plane (the contributions of the other planes cancel each other). The total electric field 
acting on these ions is then EiocP ~ E' + (k/2£o£oo) = {Eu + Eint)/2. Even this simplified picture of the electric 
fields acting on the ions allows us to explain the spectral weight anomalies. The spectacular reduction of the spectral 
weight of the 320 cm~^ phonon mode in underdoped Y123 shown in Fig. 1, e.g., is due to the fact that in the frequency 
region around the phonon the two fields Ebi and Eint have opposite signs (see Ref. iQ for a discussion). 

In reality, the displacements of the ions modify the electric fields Ebi and Eint , so that the model equations become 
slightly more complicated. As an example, we show in Fig. 7 (b), how the displacement of the planar oxygens of Y123 



(0(2) and 0(3)) from their equlibrium position influences the electric field Ebi- We refer the reader to Refs. 30 
for further details of the model 



Appendix B: Approximate "tight-binding" sum rule for bilayer compounds 

In this appendix we put the conclusions of section IV B on a more rigorous basis and we present a way of estimating 
the contribution of the countereffects We shall discuss a microscopic counterpart of the model outlined in 

Appendix A. For this reason, it may be helpful for the reader to follow Fig. 7. Let us consider a model defined by the 
following hamiltonian 

H = Hin-plane + Hu + Hint , {HA) 

where Hin-piane contains the intra-planar single particle terms and interaction terms, Hu and Hint contain the 
hopping terms between the closely-spaced copper-oxygen planes and between the widely spaced ones, respectively, 

Hu^-t±U ^ ,.Cj,;,s +H.C., {B.2) 

/— 1,3,5,... ;i;s 



Hi, 



-t±, 



E 



■-1,1+1. 



H.C.. 



(B.3) 



Z=2,4,6,... ;i;s 
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In the presence of an electric field along the c-axis, the hamiltonian reads: 



i=i,3,5... ^ Za ti / i=2Afi... ^ la h J 

Here 

i,s 

where a = bl for / G {1, 3, 5, ...} and a = int for / G {2, 4, 6, ...}, are the c-axis paramagnetic current densities. Further, 

are the c-axis kinetic energies (per unit cell), N is the number of lattice sites in one copper-oxygen plane, Abi is 
the vector potential in the intra-bilayer region {Abi — iujEbi), Aint is the vector potential in the inter-bilayer region 
{Aint = iujEint)- When studying the optical response of the system, selfconsistent values of these two vector potential 
have to be used (RPA approximation). Using similar manipulations as in Ref. | ]50| , we obtain the following relation 
between the averaged total current densities, ja , and the two electric fields, Ea , 

ja = '^a,pEi3 . {B.7) 

Here 

where I = 1 for a = bl and I — 2 for a = int. Note that all intra-bilayer regions are identical and all inter-bilayer 

regions are identical, the long-wavelength limit is assumed. Finally the four conductivities <Ja,f3 [o^jP G {bl,int}) are 
given as follows: 

0-a,/3(wj = —J- . {B.9) 

l(UJ + 10} 

Here 

/oo 

d{t - t'){[f{l, t),f{l', t')])e{t - i')e*"(*-*'^ , (B.IO) 

where I — 1 for a — bl and I = 2 for a — int; the sum runs over the odd values of I' for P — bl and over the even 
values of I' for f3 = int. We are not going to discuss the general case here. Instead we concentrate on the case where 
t±int << t±bi, i-e., the case of negligible inter-bilayer kinetic energy (fc(2)). In this case we can neglect uu,int, cfint.u 
and (Jint,int and we have ju ~ cru.biBbh jint ~ 0- The current-current correlation fmictions of Eq. (B.IO) have the 
analytic properties required for the validity of the sum rules and we obtain 



abi{Lo)dL. = -^^{kil)), (B.ll) 



where (7i,i{i-u) = ai,i,bii^)- The sum rule (B.ll) yields a formula for the change of the kinetic energy upon entering the 
superconducting state analogous to Eq. (8) or Eq. (10) {a{iv) is substituted by abi{uj) and d is substituted by du). 
For vanishingly small values of the regular part of au, we obtain A(fc(l))[meV] — —AC (wfc;[cm" -i])V4i[A]. This IS 
nothing else than Eq. (3) except for the factor of 4 Igg. In a general case, Eq. (11) yields a recipe for how to treat 
the countereffects. 
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TABLE I. Values of the fitting parameters and other numerical factors used in fitting the measured c-axis infrared conduc- 
tivity of optimally doped Bi-2212. The plasma frequency and the broadening parameter of the broad Drude peak are denoted 
by D,H and 75;, respectively. The meaning of the other symbols is the same as in Ref. [|l^. The temperatures are given in K, 
the frequencies and the broadening parameters in cm~^. The values of the numerical factors a, /3, 7 [|l^, a = 2.28, /3 = 0.64, 
7 = 1.28 have been obtained in the same way as in Ref. [^] (see also Ref. [^) using the following values of the effective 
ionic charges: e^^ = 3, ej^ — 2, Cq^ = 2, Cq^ — 2, Eq = —2. The distances between the planes of a bilayer and between the 
neighbouring bilayers are du = 3.37 A and di„t = 12.03 A, respectively. 
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TABLE II. Values of the intra-bilayer plasma frequency LJbi, the Josephson coupling energy Ej and the condensation energy 
Uo for several bilayer compounds. The values of lum are taken from Ref. jl^ (underdoped Y123), Ref. |l^ (optimally doped 
Y123), Ref. (underdoped Bi-2212); the value for optimum doped Bi-2212 has been obtained as described in the text. The 
values of the condensation energies are taken from Ref. [^^ (Y123) and Ref. [^ (Bi-2212). 
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TABLE III. Values of the parameters and other numerical factors used in modelling the c-axis infrared conductivity of 
Tl-2201. The plasma frequency and the broadening parameter of (ta and as are denoted by Qa and 7a and fis and 7s, 
respectively. The parameters of the phonons are denoted by Si, uji, 71, 52, (^2, 72, S3, cos, 73. The frequencies and the 
broadening parameters are given in cm~^. The values of the numerical factors a, /?, 7 [|18|, have been obtained in the same 
way as m Ref. (see also Ref. 0) using the following values of the effective ionic charges: e^j = 3, e^^ = 2, e^^ = 2, 
— —2. The distances between the apical-oxygen planes are dA ~ 5.3 A and ds — 6.3 A, respectively. 
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FIG. 1. (a) Experimental spectra of the c-axis conductivity, a — ai+ 102, of Y123 with Tc = 53 K. Experimental data (thick 
solid lines) and fits (thin solid lines) for (b) T — 300 K and (c) T = 4K. The symbols (A), ( B) , and (C) indicate the most 
pronounced anomalies as discussed in the text. A slightly modified version of Fig. 1 from Ref. bof. 



FIG. 2. Experimental spectra of the c-axis conductivity of almost optimally doped Bi-2212 with Tc — 91 K. (a) Data 
for T = 300K and T = 100 K. (b) Data for T = 100 K and T = 10 K. The insets of (a) and (b) display the differences 
CTi (T = 100 K) - ai (T = 300 K) and cti {T = 10 K) - ai (T = 100 K) , respectively. Note that the former and the latter difference 
corresponds to a temperature change of 200 K and 90 K, respectively. The symbols (A) and (B) in (b) indicate the most 
pronounced anomalies as discussed in the text. 



FIG. 3. (a) Experimental spectra of the c-axis conductivity of optimally doped Bi-2212 — region of the spectra near the 
anomalies on an enlarged scale, (b) Fits of the spectra obtained by using the model of Ref. jTsj as described in the text, (c) 
Spectra of ai{T = 100 K) - (Ji(T = 10 K). (d) Spectra of the quantity Nn - Ns defined in the text. 

FIG. 4. Schematic representation of the spectral-weight changes upon entering the superconducting state. The increase of 
spectral weight at low frequencies below Tc is compensated by a decrease around some frequency u* . This frequency may be 
either (a) higher or (b) lower than the interband cutoff Uc- The horizontal lines below the frequency axis of (b) indicate the 
three frequency scales corresponding to the three effective c-axis kinetic energies discussed in the text. 

FIG. 5. Results of model computations of the c-axis infrared conductivity of Tl-2201. The results shown in (a) and (b) 
have been obtained using the assignment AI (the Tl-layer oxygens vibrate at the highest frequency, the apical oxygens at 
~ 400 cm~^), those shown in (c) and (d) correspond to the reverse assignment {A2). The results shown in (a) and (c) have been 
obtained considering the region A as "intra-bilayer" , those shown in (b) and (d) have been obtained considering the region A 
as "inter-bilayer" . The planar oxygen mode manifests itself as a step around 240 cm"'^ in (b) and (d). The dashed and solid 
lines correspond to the two values of the plasma frequency of region A, Qai and ^Ia2, f^Ai < ^A2, respectively, as discussed in 
the text. The dotted and dashed-dotted lines in (d) represent results for Qa ~ 1500 cm~^ and Qa ~ 600 cm^^, respectively. 
The inset of (b) displays the step in the experimental data of Ref. . The inset of (c) displays the experimental spectra of 
(A^„ — Ns)/ps taken from Ref. The arrow demonstrates the increase of the spectral weight of the phonon peak at 390 cm~^ 
in the superconducting state, as discussed in the text. 



FIG. 6. Crystal structure of Tl2Ba2Cu06. 

FIG. 7. (a) Schematic representation of the model, (b) The average electric field Eti between the charged planes correspond- 
ing to the planar oxygens (thick horizontal lines) displaced from their equilibrium positions (thin horizontal lines) possesses a 
contribution AEbi due to the charged Y-plane (dashed line), AEm = [((di,i/2) + a)EY + {{dbi/2) — a)[—EY)]/ dti- Here Ey is 
the field due to the Y-plane in the upper region of the figure. A slightly modified version of Fig. 2 from Ref. pol] . 
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